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Abstract 

We present an analytical framework for modeling a priority-based load balancing scheme in cellular 
networks based on a new algorithm called direct retry with truncated offloading channel resource pool 
(DRk)- The model differs in many respects from previous works on load balancing. Foremost, it incor- 
porates the call admission process, through random access. In specific, the proposed model implements 
the Physical Random Access Channel used in 3GPP network standards. Furthermore, the proposed 

m . 

I/"") . model allows the differentiation of users based on their priorities. The quantitative results illustrate that, 
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for example, cellular network operators can control the manner in which traffic is offloaded between 

neighboring cells by simply manipulating the length of the random access phase. Our analysis allow to 

quantitatively determine the blocking probability individual users will experience given a specific length 

of random access phase. Furthermore, we observe that the improvement in blocking probability per 
• i-H , 

shared channel for load balanced users using DR/< is maximized at an intermediate number of shared 

H ' 

channels, as opposed to the maximum number of these shared resources. This occurs because a balance 
is achieved between the number of users requesting connections and those that are already admitted to 
the network. 
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I. Introduction 

Given the rapid current and expected growth in HSPA and LTE-based networks and in the number 
of mobile devices that use such networks to download data-intensive, multimedia-rich content, the need 
for QoS -enabled connection management is vital. However, the non-uniform distribution of users and 
consequent imbalance in usage of radio resources leads to an existence of local areas of under- and 
over-utilization of these resources in the network. This phenomenon results in challenging network 
management issues. Load balancing is an important technique that attempts to solve such issues, and 
occurs when a centralized network controller intelligently distributes connections from highly congested 
cells to neighboring cells which are less occupied. This allows for an increase in network subscriber 
satisfaction, as more subscribers meet their QoS requirements. Furthermore, it allows for an increase in 
overall channel utilization. 

Load balancing is popular among the academic community and has been explored for many years 
as described, for example, in earlier works such as 0, [3J, as well as more recently in (H-Q, and 
many other papers. For simplicity, in the previous treatment of load balancing, it has been assumed that 
the connection admission process can be neglected. In cellular networks, however, each new connection 
needs to first send a request to the serving base station (BS) through some predefined control channel. 
In 3GPP standards, this control channel is the Physical Random Access Channel [8, Sec. 2.4.4.4], (9j 
Sec. II]. The success of a connection request by a user is dependent on multiple factors including the 
number of requesting users; the pairwise channel quality between the user and the serving base station 
(measured, for example, in BER or outage probability); and the actual control channel access technique 
itself |[T0| - lfl2ll . Until now the impact of random access overhead on load balancing performance is not 
well understood. More specifically, the quantitative relationship between the random access phase length 
and the user blocking probability in a load balancing-enabled cellular system is unknown. 

The present paper aims to provide a framework to address this gap. It provides a set of associated 
results that describe the benefits of load balancing in cellular networks in high detail, using realistic 
scenarios and network configurations. At the highest level, the approach is based on a Markov model in 
which the state at each point in time provides a full description of user connectivity, and the transitions are 
determined by a range of factors described below. While many previous studies have also applied Markov 
methods to model network traffic, the approach described here is differentiated from previous work in 
multiple respects. First, as noted, it incorporates the call admission process, i.e. the process of requesting 
a connection before the actual assignment of a BS resource. Second, specific consideration is given to 
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channel quality information, which in turn has a direct effect on the number of successful connections 
and terminations thereby having an impact on load balancing performance metrics. Third, it incorporates 
priority differentiation in the load balancing process. We use the prioritization approach, where users 
that do not have access to multiple base stations get a higher priority in using available cell resources 
than users that do. More generally, the model allows substantial flexibility with respect to traffic types, 
admission control, densities of user population, and other parameters, but is also specifically constructed 
with the function of considering optimal network policy decisions with respect to load balancing. 

The rest of the paper is organized as follows. The system model is introduced in Section QI1 while the 
analytical model is introduced in Section JII] The numerical results are presented in Section [IV] Lastly, 
the paper is concluded in Section [Vj 

II. System Model 

A. Channel Structure 

We consider a cellular system where two base stations (BSs) are positioned such that they create a 
region of overlap in coverage. Naturally, while cellular systems typically have far more than two BSs, 
a reduction to a two-BS system for analytical reasons enables a tractable analytical framework while 
still allowing exploration of a large number of microscopic parameters to use in optimizing network 
performance. As our model is of microscopic nature, the extension to a multi-BS system is beyond the 
scope of this paper. We strongly emphasize that consideration of load balancing in the context of a two- 
BS system only has been used extensively and successfully in previous treatments, e.g. (U, Q, lfT3l , 

m. 

Cell 1 has M\ available basic bandwidth units, as referred to in lfl31 or more commonly referred 
to as channels, and cell 2 has M2 available channels. The throughput of every channel in each cell is 
the same and equal to //bits/second. We assume that channels are mutually orthogonal and that there 
is no interference in the set of channels belonging to cells 1 and 2. Each BS emits a signal using 
omnidirectional antennas and we assume a circular contour signal coverage model, in which full signal 
strength is received within a certain radius of the BS, and no signal is available beyond that radius, as 
used in e.g. @, 0, lfl6l 

B. Node Placement 

Each mobile terminal, referred to as user equipment (UE), remains at fixed positions following an 
initial UE placement process, with each UE located in one of three non-overlapping regions as shown 
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in Fig. Q] N\ UEs are in group 1 and have access only to one BS, N2 UEs are in group 2 and have 
access only to the second BS, and N3 UEs are in group 3 and can potentially access either BS. Such a 
non-homogenous UE placement has been considered, for example in Ifl4l . which allows for a tractable 
analysis of the considered system and considers all important groups participating in the load balancing 
process. The non-homogenous case of UE distribution is the most widespread because UEs are generally 
distributed non-uniformly over a cellular area. UEs in group 3 are in the region of overlap in coverage 
between the two cells, also known as the Traffic Transferable Region (TTR) (H. Since only one UE can 
occupy a channel at a time, a maximum of Mi + M2 UEs can be connected to both cells in the system 
at any given time. We assume that UEs from groups 1 and 3 are initially registered to cell 1 (serving as 
the overloaded cell), while UEs in group 2 are initially registered to cell 2. 

C. Signal Transmission Model 

Time is slotted and the minimum time unit is a frame length of t s seconds. Connections and channel 
conditions are assumed to remain constant for the duration of a frame, though they will in general vary 
from frame to frame. Moreover, we assume that the connection and termination processes for UEs are 
directly dependent on the channel states experienced between each group of UEs and the BSs they are 
connected to. We also assume that channel states are binary and independent from slot to slot, much as 
occurs in ifTTI . and that all of the UEs in each group experience the same channel quality to a given 

(i) 

BS. Therefore, in any given time slot a UE is either experiencing a good state with probability w x ' y , (x 
denotes the particular pair-wise connection between group x of UEs and associated BS y, and i G {d, u} 

(i) 

denotes the downlink and uplink, respectively), or a bad state, with probability 1 — w x ' y . The value of 

(i) 

Wx'y is dependent on the distance between UEs in group x and BS y, which is denoted as d XtV . In our 
analysis we use the distance as an input to a combined path loss and shadowing propagation model. 

D. Prioritization 

Because of the strict boundaries between groups of UEs, we assign priorities on a per group basis. A 
single, higher priority is given to all UEs in groups 1 and 2 because there is no flexibility to reassign 
them to a different BS; group 3 UEs can potentially be reassigned and thus given a lower priority. A 
very similar priority model has been used in other treatments of load balancing. For example, in [13], 
lfl"8l . newly arriving connections in the non-TTR are given first priority to acquire channels from their 
serving BSs, while the connections from the TTR are assigned to the remaining channels. Our model 
allows for the implementation of a wide range of scenarios that require such traffic prioritization. One 



Fig. 1. The system model considered in the analytical evaluation. The two large circles represent two adjacent cells, cell 1 and 
cell 2, containing 3 groups of UEs. Each UE is indicated by a small circle, and the triangles represent the BSs at the center of 
each cell. UEs in group 1 have access only to the BS of cell 1, UEs in group 2 have access only to the BS of cell 2, and UEs in 
group 3 have access to both BSs. The region of overlap between the two cells is referred to as the Traffic Transferable Region. 

potential application is for the modeling of networks where load balancing traffic originating from UEs 
in the TTR has lower priority than non-balanced connections due to several factors including a lower 
average channel quality fl9l , [20], QoS requirements 11211 . non-uniform spatial distribution of traffic 
classes lfl4l . or cell dwell time. Furthermore, it allows for modeling integrated hybrid cellular/WLAN/Ad 
Hoc networks as discussed in @, 11221 . where non-cellular terminals in the TTR have a lower priority than 
cellular UEs, and hierarchical cellular systems [23], where members of different tiers have independent 
priorities. Finally, it enables the modeling of femtocell traffic prioritization, where UEs in groups 1 and 
3 are those in the Closed Subscriber Group (CSG) E4l . while UEs in the TTR are neighboring UEs 
outside of the CSG. 

E. Random Access 

In the connection process a UE first attempts to connect to the BS it is initially registered to by 
requesting a connection through a random access channel. We assume a frequency division duplex 
transmission mode, where control and data traffic are transmitted and received simultaneously. Each 
UE generates a connection request with probability p x . A connection is requested randomly in one of 
L x < t s non-overlapping, time slotted control resources, unique to group x G {1, 2, 3} of UEs. In other 
words, each group has a unique set of sub slots within a frame during which UEs may, but are not 
required to, request a connection. The random access phase length is equal to slot length t s . Collisions 
between connection requests from UEs in the same group are possible. 
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The random access procedure considered in this work shares features of the 3GPP-based cellular 
network standards, which use the Physical Random Access Channel (PRACH), mapped on a one-to-one 
basis to the logical random access channel (RACH). RACH uses the S-ALOHA protocol and, in relation 
to the priorities assumed in this paper, allows the prioritization of connection requests based on Active 
Service Classes (ASC) |[T0l which are unique to each UE, and can be adapted by the 3GPP-MAC layer 
once the UE are in connected mode JU Sec. 2.4.2.6]. The BS advertises itself to the UEs within range 
through the broadcast channel using signatures (3GPP release 99, e.g. UMTS), subcarriers (3GPP release 
8, e.g. LTE), or time slots, which each ASC can in-turn use for connection requests on RACH. The 
adaptation of the ASC is performed in the time intervals predefined by the operator. For the purpose of 
our paper we assume that the BSs collectively, through the Radio Network Controller, map the received 
signal from every registered UE to an associated ASC. 

We assume a zero-persistence protocol, i.e. a collision during a connection request implies that con- 
nections are lost, and also UEs do not retry to generate another dependent connection. Due to this 
assumption a power ramping process, i.e. feedback from the UE to the BS on an unsuccessful connection 
request ifTOl Sec. ITB], is redundant. To isolate the impact of each group of UEs on collision rates, we 
assume mutually exclusive RACH resources assigned to each ASC Qj, Fig. 4]. Analysis of PRACH 
performance in isolation can be found in |9j, ifTTTl . 

F. Data Transfer 

A connection request is granted during the connection arrangement process if a good channel state 
occurs between the UE and its associated BS at the moment of the request, and if no collisions occur 
between multiple requests from different UEs. Once a connection is established, the BS randomly selects 
a channel and assigns the connected UE to it. The UE then begins to receive downlink data. UEs terminate 
a connection with probability q, where 1/q = r p /(Ht s ) is the average connection transfer size and r p 
is the average packet size given in bits. We assume that the transfer size is at least one frame long. A 
connection terminates either when a transmission completes, or when the channel is in a bad state during 
transmission. 

G. Load Balancing Process 

In the case of a UE in group 3, if a connection request is successful and there are no resources available 
in cell 1, we assume that the radio network controller performs load balancing by transferring the call 
from cell 1 to cell 2. To avoid overloading cell 2 and protecting UEs that are already registered to it, UEs 
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in group 3 can access a maximum of K channels from cell 2, where < K < M 2 . UEs in group 3 have 
access to an additional K shared channels, therefore they have access to a total of M\ + K channels. 

Using the nomenclature of (4j Sec. 2 and 3] this load balancing scheme belongs to a class of direct 
load balancing schemes. It is closest in operation to direct retry 0. Since we allow at most K available 
channels to offload traffic from cell 1 to cell 2 (as in simple borrowing scheme \4, Sec. 2]) we denote 
our scheme as direct retry with truncated offloading channel resource pool (abbreviated as DRx). With 
K = M 2 our scheme reduces to classical direct retry, while for K = it reduces to a system in which 
no load balancing takes place. 

In our model we do not use a take-back process. That is, once connections from group 3 are offloaded 
onto cell 2, they remain connected to cell 2 during the transmission despite whether or not resources 
have been freed in cell 1. In |5] the authors remark that the take-back process, although more fair to cell 
2 because it minimizes blocking at cell 2, is not always advantageous to the network due to the high 
signaling load that accompanies it. Additionally, as in |@), lfT6l . we do not use queuing, so there is no 
consideration of a call give-up process (5). Moreover, we do not allow preemption of connections from 
the TTR by connections that have access to channels only from their respective BSs. 

Note that all variables introduced in this section, as well as all other variables used in this paper are 
summarized in Table U 

III. Analytical Model 

Let {A, Y^ 1 ' , W 2 ), C} denote a state of a Markov system, where A denotes the number of resources 
used by group 1 UEs, Y^ 1 ' and denote the number of resources used by group 3 UEs associated 
with cell 1 and 2, respectively, and C denotes the number of resources used by group 2 UEs. Then the 
steady state probabilities can be denoted as n a ,b,c,d — Pr(^4 = a,Y^ = b, Y^> = c,C = d). Note that 
a+b+c+d< mm{N 1 + N 2 + N 3 ,M 1 + M 2 }, a + b < min{iVi + N 3 , Mi}, c + d < min{N 2 + N 3 , M 2 }, 
and b + c < minjA^Mi + K}. These conditions govern what states are possible in the transition 
probability matrix. 

We define the state transition probability as 

(a t ,b t ,c t ,d t ) in /1 _ n _ l y(2) _ n — A 

\At-i = at-iM-i = h t-iM 2 \ = ct-i.Q-i = dt-i), (1) 

where subscripts t and t — 1 denote the current and the previous time slots, respectively. This allows 
for computation of the transition probability matrix required to obtain n a 6 c d, which is in-turn used to 
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TABLE I 

Summary of the Variables Used in the Paper: Functions (Top) and Variables (Bottom) 



Variable 


Description 


Units 


{A}, {Y (1) ,r (2) j, {Cj 

L J ' L ' J ' J 


Number of channels used: group 1, 3 (connected to BS 1 and BS 2, resp.) and 2 




7T r, 7i rl 


Steady state probability of the considered Markov chain 




(ot,6 t ,c t ,d t ) 
a t-l. 6 t-l. c t-l> d t-l 


Transition probability of the considered Markov chain 




('»qS'! (-)fl(-) (•)£(•) 

"(•)' (•)' MO 


Functions suDportine T X a tM,c t ,dt) 




-P w (■) 


Distribution of 7 at distance d XjV for i = {d,u} 




(*) 


Outage nrnhnhilitv bptwppn trrniin t and RS ?/ in ?' — -f*7 

LlLtlii^. IJ1 *J L>CLL>1±1 L V gLUUU J. ClllU _UO U 111 t \ 1 Lt f 




O(-) 


Q function and indication function, resp. 






Blocking and collision probability for group x, resp. 




qU) rp(j) 


PrrVnd ni 1 1 1"\7 rvr ru^ ti »- q 1 1 ti <t nnd fprtninntinfr n cfinnprlinnc <r i ^ rp r\ r> acti'vc* licf^vc rpcn 
rlUUdUlilly yJL gCHClaLlllg, ilUU LCIllllllaLlllg, J L-UllllC^UUlls glVCll i aULlVC UaCls, ICap. 




oU) 
l J k 


Prnhabilitv nf i" snrrps^fnl rnnnppHnns nut nf ^ rpniipst*; nn R Af^l-T 

L lUUttUlll IV V / 1 J 3 Ul_\yLj l3± Lll \^W1111V^V^L1»J11l3 KJ LI L V / 1 f\i I^IJ LIAjJ V /1 1 l\ /VV 1 1 




,U) 

i,x,y 


^liinnorf in <t fnnption fo rnmnnfp ofj) 

O UpjJUl Llllg 1 LlllL. LlUll L(J ^(JllljJLlLC i_7^ 




L x , L 


Number of access slots assigned to group x, and all groups, resp. 






TnHp"YPQ Hpnntincr nnlinV anH HownlinV rpt;n 




y 11 ?) i i 1c 

A, U/, [V , (j, /, rv 


Snnnnrtin (t variflhlps 

iJ UUUUl Llllg V LlLlCHJl\sJ 




(X, &, C, 


annpl c or'pnni pH * crrnnn 1 crrr\nn ^ ( pnnnpptprl f o R^l 1 anH RQ rpsn ^ ornun 
\^.lld.llllCla (JCUU.U1CQ. ULIU 1 , iilULlU J ^UUllllCC LCU. L(J Do 1 ill 1U 13 O Z,, ICaU.J, gl UUU Z. 




fc, t — 1 


Current and previous time slots indexes, resp. 






Frame length 


seconds 


(*) 

7q 


oignai recepuon inresnoiu ior z — -^ix, tjj- 


Qd 


7 


D pr>pll rp H c 1 Crn Q 1 1 

JvCCClVCU ai^Ilal 1CVC1 


rlR 
UD 




Ptittilrwc PYnnnpnl onH cfmrln\s7i n cr \?^\r^ dtipp rpsn 
rallllUsa CApUilClll clllU allaQUWlllg Vallall^C, ICap. 






1? cic i \ 7(± r\ c i rrn a 1 at o" — J rl o i \. "Frvr T TTh rrrnnn ~r q ti H RQ n 
JxCCClVCU ol^Ilal ill i — 1*^5 f A ^ Uli glULip X allli Do y 






Prnhfibilitv of pnnnpptinn rpniipst (rpnprfltinn bv TIP Qrnim t and all (rrnnns rpsn 

X lk_Jl_Jcll_JlllLy V / 1 ^Wllll^w L1W11 ILUUUk^L CILIWII LJ V V_J 1 _. ill ^ LI VJ Jb , CL11L1 CL11 clUUL/iJ^ IWtjLJ. 




d d (i) 


Normal and reference distance for i = {w, t/j- used in pathloss calculations, resp. 


m 




Individual rhannpl thrniiffhrnit 

1 1 1L11 VluUul CllClllll^l LlllWLI.gllL^LI.L 


b/s 




Prupfinn nf pViannpl^ in R^l / 




JJ 


f^Vitinnpl ntili7titi on 




w w 


Unit-less constant used in pathloss calculation for i — {u,d} 




p(0 


Transmitted power by BS (i = d) and UE (i = u) 


dB 


m h , at* 


Number of channels used by BS y and number of UEs in group x, resp. 






Average packet length, given channel quality transmitted, from BS y to UE group x 


frames 




Average packet length in perfect channel conditions 


frames, bits 




Required probability of blocking and collision, respectively, for group x 




Pa, C/c, Pd, Pe, io, 


Variables and functions supporting 




*7 


Variable supporting 
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compute the performance metrics of the load balancing-enabled cellular system. In the subsequent sections 
we describe the process of deriving the transition probability '^''^ 1 d t 1 - We begin by explaining 
the computation process for the channel quality, and then focus on the derivation of the functions that 
support (fl}. 



A. Derivation of Channel Quality 

In the downlink, the probability of a UE belonging to group x 6 {1, 2, 3} and receiving a good signal 
from BS y € {1, 2}, is defined as 



w%l= i - w!",; £ -:; ,! ) ■ .1 / /v (-.</, ,,)<!-. ^ 

where 7^ is the signal reception threshold for the downlink, expressed as the minimum required received 
power of the received signal J^y, and p s w (7, d x>y ) is the distribution of the signal 7 received at group 
x, which is at a distance of d XtV from BS y. As an example, we consider an environment with path loss 
and shadowing, where w x d } y is expressed in 11251 Eq. 2.25] as 

jd) _ p (d) _ lologio w (d) + 10< yi ogl A 



A 



1 - Pr(45 < 7 f°) = 1- / 9 p s wh,d x , y )dj, 

Jo 



% - n - wio gl0 w v y + iuo iog 10 ^ 



wffi = Q I ^ I , 0) 



where is a unit-less constant dependent on both the antenna characteristics and an average channel 
attenuation, given as |[25l Eq. 2.41], is the BS transmitted power (which is assumed to be the same 
for both base stations), d XiV is the distance of the UE in group x located furthest from BS y; d^ \ is the 
reference distance for the BS antenna far-field; c<j, is the log-normal shadowing variance given in dB; 5 
is the path loss exponent; and the Q function is defined in the usual manner as 

Q(z) 4 / e - u2 ' 2 du. (4) 

Jz V27T 



(u) 

In the uplink, the probability that a good signal is received by BS y from a UE in group x is w X} y, and 
is expressed in the same manner as equations Q and (O, by replacing all variables having superscript 
(d) with (u), where pj u ' denotes the UE transmission power; denotes the constant for the UE 

antenna; d^l is the reference distance for the UE antenna far-field; and 7,^ is the signal reception 
threshold for the uplink. The downlink and uplink channel quality information governs the success rate 
of the connection admission process, as well as the duration of a downlink transmission. 
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B. Derivation of Connection Arrangement Probability 

An important feature of the model is the consideration of connection admission in the load-balancing 
process. This process is a function of the total number of UEs, the number of UEs receiving data from 
their respective serving BSs, the pairwise channel quality between the UEs and its serving BSs, and the 
underlying random access algorithm. The probability that j new connections have successfully requested 
downlink data given %t-\ G {at-i, &t-i, ct—\, dt-i} currently active connections from group x, with a 
random access channel consisting of L x time slots is 



c(i) A 



T (it) 7 > 



1 _ ,,,(«) + J 1 ) j - n (5) 

w x,y T" T,; x , J — U, W 



X >V ^ i t -i,x,yi 

0, otherwise 



where 



N a -i t 



#U = «S E~ ^ x )rfd - fc)^- 1 "^; w 

Pa; is the probability of a connection request by an individual UE in group x; and BY 1 is the probability 

that among k UEs requesting a connection, j were successful in obtaining a resource. Note that the 

(i) 

reference to x,y in SY is omitted due to space constraints, keeping in mind that for at,at-i x = 1, 
y = 1, for &t, 6 t _i x = 3, y = 1, for q, q_i x = 3, y = 2, and for c^, c^_i x = 2, y = 2. 

For consistency with cellular networks such as 3GPP, we consider a PRACH-like control channel, for 
which can be described in the manner of ||26l Eq. (3)] 



min(k,L m ) 

(m — j)\(L x — m)\(k — m)\n\ 



B U)_ V (-l) m+ HL x -m) k - m k\ 



(7) 

Note that depending on the assumption of how collisions are resolved, different definitions of BY in 
equation ([7]) can be applied when calculating the connection arrangement probability according to ©. 



C. Connection Termination Probability 

Once a UE successfully requests a connection from the serving BS, a downlink transmission is started 
provided that at least one free channel is available for the UE. The connection terminates when the BS 
finishes transmitting data to the UE or when the downlink signal received by the UE is in outage. The 
probability that j connections from i t -i active connections at group x terminate is 
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where l x>y = 1 — w^y + w^ y q is the inverse of the average packet length, accounting for truncation 



of some packets due to a bad channel quality. Again, the indices x, y have been omitted for notational 
simplicity in the syn 
Section InEE holds. 



simplicity in the symbol 1^ _ , assuming that the same relationship between x,y and j,it-i as given in 



D. Derivation of the State Transition Probability 

Using the definitions of the arrangement and termination probabilities, expressed in ((5]) and ([8]), 
respectively, we can finally introduce the state transition probabilities for the complete model. The 
transition probability is constructed using the termination and arrangement probability definitions and 
the respective relationship between the variables of these two definitions (which are dependent on the 
start and end states of the transition). Due to the complexity of the derivation we begin with a highly 
simplified example. 

1 ) State Transition Probabilities for a Single UE Group: To facilitate understanding the derivation of 
the complete state transition probabilities, we first consider a network in which no load balancing occurs 
and all of the UEs are in group 1, such that Ni > and N2 = N% = 0. The state of the Markov chain 

(a f ,0,0,0) 
'ot-ijO.O, 1 

EZoTt^-^sZL <H-l > a t , at / M u 

T,T=^l 1 Sit a r at - 1 \ ot-i < at, at ± M U 



then simplifies to {^4, 0, 0, 0} and the transition probability becomes r^QQQ, where 



(ot,0,0,0) 
r a f _ l5 0,0,0 



(9) 

Eat rp(i) q{i+a t —a t -i) 
i=0 J a t _i»Ja t _i 

1 v^a* rp(i) v^JVi o(i+j-at-i) ^ _ M 

+ l^i=0 1 a t -i l^j=M 1 b a t -i , Ot-l < a t , a t — Mi. 

In (O the case a t -\ > at, at 7^ M\ denotes the transition from a higher to a lower channel occupancy, 
subject to the constraint that the number of channels occupied in the end state can not exceed the total 
BS capacity. The number of terminating UEs is set to compensate for the UEs that generate successful 
connections. The case at-i < at, a t 7^ M\ denotes the transition from a lower to a higher channel 
occupancy (given, again, that the number of occupied channels is less than the total BS capacity). In 
this case UEs from group 1 need to generate exactly as many connections as given by the end state, not 
forgetting to generate connections in order to compensate for the total number of terminations. Lastly, 
for the case of at-\ < at, a t = M\ the end state equals the total channel capacity. The first term in 
the definition of this transition probability includes exactly the number of connections needed to reach 
the end state, once again compensating for terminations. The second term accounts for all successful 
connections generated that exceed those needed to reach the end state, which will not be admitted. 
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2 ) General Solution for the State Transition Probabilities: As expected, the current state of the UEs in 
each group has a strong impact on the allowable future states in all groups. For example, newly generated 
connections from UEs in group 2 are initially assigned to cell 1 and can only be offloaded to cell 2 once 
all channels on cell 1 are occupied. Therefore, in the process of constructing the complete set of transition 
probabilities, not only is it necessary to enumerate all conditions that govern these possible transitions, 
as in (O, but it is also important to consider the relationship between the possible combinations of 
termination and arrangement probabilities. 

Before presenting the general solution, we introduce supporting functions that simplify the description 
of the transition probabilities. First, because the system is composed of three groups, where groups 1 
and 2 have the same level of priority, we define a function that governs the transition probabilities for 
these groups, i.e. 



(k+it-i-it) q(i) 



(fc) Q (k+i t -i-it) 



(it,z) . 



k=0 1 i t . 1 D i t . 1 



Eit rp(k) q(k+i t —it-i) 
fc=0 1 i t -iit-i 



h + jt < Mi, i t < i t -i', 
it+jt < Mi, i t > i t -i or 

it + j t = Mi, i t > it-iJt > jt-i\ 



(10) 



jt-l—Z+lt 



it-i 2—in=i t +l %t-\ 



(k+n— it-i) . 



; i t + jt = Mi , i t > i t -i , jt < jt-i- 



where I, 



U) 



r(i) 



1 when i > j and If = 0, otherwise. Variables i x and j x are supporting parameters that 



will be replaced by respective variables of £0), once we derive general formulas for transition probabilities. 
The ranges of i x and j x will be defined in the respective transition probabilities shown later. Note that 
(fTOl) resembles (|9]), except for the introduction of the indicator function If\ For the remaining groups, we 
define the following supporting functions which denote the possible termination probabilities for group 
3 



(n)n( b t,ct,k,l) ^ 
6 t _i,c t _i ^ 



rp(k) rp(l) 


n 


= 1, 


rp(k+b t -i— b t )rp{l+Ct-x-C t ) 


n 


= 2, 


rpik+bt-!-^) rp(l) 


n 


= 3, 


rp(k) rp(l+Ct-l-C t ) 


n 


= 4. 



(11) 



Note that the termination probabilities for group 3 in (fTTT) are composed of individual termination function 
as given in ([8]). One reason for this is that different channel qualities are experienced by UEs in the TTR, 
resulting in unequal termination probabilities, depending on whether these UEs are connected to BS 1 
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or BS 2. Lastly, we define 



(n) Jb t ,c t ,k,l) A_ 
?6 t _i.c t -i ~~ 



(12) 



^(k+l+bt+Ct-bt-i-Ct-!) _ 1 

Q (fe+/ + C t -C f _!) Q 

ot-x+Ct—i 
5 (JH+J-*.-0 n = 4 . 

Ot-l+Ct— 1 

which denotes the possible arrangement probabilities for group 3. Note that the variables k, I and z 
in ([Tol l. (fTTT > and (fT2l are the enumerators. Given the above, we can identify two major states of the 
system as follows: when all channels are occupied on both cells (called an edge state and having the 
same meaning as the last condition in (©); and the remaining states. We start by describing the edge 
state conditions. 

3) Edge state: Here we list the following sub-cases. For at + bt = Ml, q + d t = M 2 , h > &t-i, 
ct > ct-i or at + b t = Mi, c t + d t < M 2 , b t > bt-i, c t > c t -i, c t = K we have 

b t -i c f _i 

(a t ,b t ,c t ,d t ) _ \ ^ \ ^ / (l)a(bt,ct,k,l)(i)Jb t ,c t ,k,l) (l) „(b t ,Ct,k,l) 

a t -i,bt-i,ct-i,dt-i / j / 4 \ h-i,ct-i >6t_i,Ct_i b t _i,c t _i 
fc=0 Z=0 



(l)^(6 t ,ct,fc,r)^ ^ A , ,^( a " fc ) ,„„ „ ,„(*,0 
r=/+l 

Equation (fT3l holds when the number of connections from group 3 UEs to both BS 1 and BS 2 increases. 
The first term in the brackets enumerates all the possible cases of terminations and generations in group 3 
given a certain starting state. The second term in the brackets accounts for the edge case. This condition 
is similar in nature to the third case in ©. Lastly, the remaining ^afrl terms account for the possible 
transitions in group 1 and 2. The indicator function used in the last condition of (fTOl ) is a function of 
the termination and connection enumerators in group 3. That is, depending on how many connections 
are admitted in BS 1 and BS 2 in a previous frame, a certain number of UEs from group 1 and 2 that 
request connections will not be admitted. 

For at + bt = Mi, ct + a\ = M 2 , b t < b t -i,c t < c t -\ or a t + b t = M\, c t + d t = M 2 , b t = 
bt-l,<h < c t -i or a t + b t = Mi, c t + d t = M 2 , b t < bt-l,<h = Q-i or a t + b t = Mi, c t + d t < M 2 , 
b t < bt-i,ct = ct-i = ifwe have 

b t c t 

(a t ,bt,ct,dt) _ ' 

at-i,6t-i,ct_i,dt-i ~~ 

k=0 1=0 

X 

r=l+l 



(a t ,bt,c t ,d t ) _ ST~^ ST~^ ( (2)n( b t,c t ,k,l)(2)Abt,c t ,k,l) (2) n{b t ,c t ,k ,1) 

Ot-i,6t-i,ct_i,dt-i ~ / j / j V 6 t _i,c t _i ^6 t _i,ct_i * 6 t _i,ct_i 

E(2) e (fet,c t ,fc,r)\ (a t ,fc+6 t _i-6 t ) (d t ,/+c t _i-c t ) n 



14 



In this case the number of connections from group 3 to both BS 1 and BS 2 decreases, or the number 
of connections in any one of the BSs remains the same, while the other decreases. The construction 
of the transition probability is the same as in ([TBI , respectively replacing ^9^1 with ^9^}. Note that 
the definition of ^9^1 in ([Til defines the number of freed connections at BS 1 and BS 2 because of 
terminations of group 3 UEs. Since the number of connections have to be maintained at full cell occupancy 
for cell 1 and cell 2, the respective functions for BS 1 and BS 2, are used to compensate for the 

possible number of terminations at each cell due to group 3 UEs in order to maintain full system-wide 
occupancy, i.e. to remain at the edge state. 

For a t +b t = Mi, c t +d t = M 2 , b t < b t -i,c t > ct-x or a t +b t = Mi, c t +d t < M 2 , b t < h-i, °t > <h-x, 
c t = K 

b t c t -i 

(a t ,b t ,ct,dt) = ST^Sr^ f{3)n0>t,ct,k,l)(3,)jbt,ct,k,l) (3)a(h, c t ,k,l) 

k=0 1=0 

N 3 

v , (3) Ah,c t ,k,r)\ (a t ,k+b t -i-b t ) (d t ,l) , AC -s 

X ^ f -iA-i j (M^K^,^ (M 2 ,N 2 ) a d t _ 1 \c t - 1 ' ( 15 ) 

r=/+l 

This case describes the situation where the number of connections from group 3 UEs to BS 1 strictly 
decreases, while those from group 3 UEs to BS 2 strictly increases. The transition probability represented 
in (TT3l can account for this by replacing ^9^1 with ( 3 ) 9^}. The definition of ^6^1 from (fTTT> describes the 
case when the number of terminations at BS 1 from group 3 UEs account for the decrease in the number 
of connections, while the number of terminations at BS 2 account only for any additional number of 
generations. The respective function (.)a£j for BS 1 and BS 2, again, must compensate for the changes 
in connections from group 3 to BS 1 and BS 2 in order to maintain full-cell occupancy. Lastly, for 
a t + b t = Mi, ct + dt = M 2 , b t > b t -i,c t < c t -i 

bt-i c t 

(at,bt,c t ,dt) _ ST^ ( (4) fl(6 t ,ct,fc,0(4) Ab t ,c t ,k,l) (4) n{b t ,ct,k,l) 

o t _i,6 t _i,Ct_i,d t _i / j 7 / [ 6 t _i,c t _i >6 t _i,c t _i ' 6 t _i,c t _i 
fc=0 1=0 

N 3 

X ^ftt-LCt-! J (M 1 ,W 1 )«a f _ 1 ,6 f _ 1 (M 2 ,JV 2 )« di _ 1 , Ci _ 1 • (16) 

r=/+l 

The above case is the opposite of that in (fl~5l ). Here, the number of connections from group 3 to BS 
1 strictly increases, while those from group 3 to BS 2 strictly decreases. Again, respective expressions 
for Wflj;] in (TO]) need to be replaced by ^9^. The explanation for ^9^ and given in ([111 is 

equivalent to the explanation for (TT5T ). 

Non-Edge state: The second major group of cases refers to the situation in which the number of 
connections at BS 1 or BS 2 is less than or equal to the maximum capacity. This obviously involves 
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more cases to consider than those explained in Section IIII-D3I We start by denoting conditions under 
which a transition from one state to another is not possible. That is, for at + bt < Mi, ct + dt < Mi, 

b t > b t -i, c t > ct-i or a t + b t < M x , c t + d t < M 2 , b t < b t -i, c t > c t -i 

(a t ,b t ,c t ,dt) „ (M\ 

r a t - 1 ,b t -i,ct-ud t -i ~ U - KLI) 

For a t + b t < M x , c t + d t < M 2 , b t > b t -i, c t = c t -\ 

b t -i 

(a t ,b t ,ct,dt) _ ST^ (1) a(bt, c t , k, 0)(1) J b t ,c t ,k,0) (a t ,k) J dt ' ^ (\ 9.\ 

r at-ubt-i,c t -i,dt-i ~ 1^ °bt-i,ct-i Sbt-uct-i (Afi.JVi) a tt t _ 1 ,6 t _ 1 (Mi,iV a )«d t _ ll c t _ 1 - ^> 
k=0 

The above case is partially equivalent to (fT3T ) and considers the situation where the number of connections 
of group 3 UEs connected to BS 1 increases and those to BS 2 stay the same. Also, the number of new 
connections at BS 1 is less than its maximum capacity. Since this is not an edge case for the system, 
an additional third summation term is not needed as seen in (TT3T ). Note that this condition only contains 
one summation because the number of terminations from UEs connected to BS 2 cannot exceed the 
resultant connection state. This is because if they do exceed the desired number of terminations, UEs 
that generate connections to compensate for additional terminations will instead choose to connect to BS 
1 (the BS they are registered to) thereby changing the resultant connection state. Now, for at + bt = M\, 
ct + dt < M 2 , b t > b t -i, c t > ct-i 

bt-i c t -i 

(at,bt,ct,dt) _ \^ST^ (l)a(bt,Ct,k,l)(i)Jbt,c t ,k,l) AfH,k) Ad t ,l) nQ . 

r a t -i,bt-i,ct-i,dt-i -1^1^ \-i,c t -i Zbt-i,ct-i (M^OViA-i^A)^.,,^,- (19) 
k=0 1=0 

This case is an extension of the case described in (fT8l ). However, full-cell occupancy now occurs at BS 
1, i.e. all channels of BS 1 are occupied after the transition, and BS 2 operates at less than its maximum 
capacity. An additional summation is used as compared to (fT8l ) because full-cell occupancy on BS 1 
allows for terminations to occur on BS 2 from group 3 UEs without changing their resultant connection 
number. Now, for a t + b t < Mi, c t + d t < M 2 , c t < c t -i, b t < b t -i, or a t + b t = M\, c t + d t < M 2 , 
b t = b t -i, c t < ct-i or a t + b t = Mi, c t + d t < M 2 , b t < b t -i, c t = c t -i 



(a t ,b t ,Ct,dt) ST~^ ( 2 ) Q (bt,c t ,k,0)( 2 ) Jbt,c t ,k,0) (a t ,k+bt-i-b t ) (dt,c t -i—c t ) /om 

r a t -i,b t -i,ct-i,dt-i =2^ d bt-i,ct-i ^-i.Ct-1 {M u Ni)U atiMi (M 2 ,N 2 ) Ct dt-i,Ct-i • ( 2U ) 

k=0 

The case described by d20l ) is a direct extension of (fT4l ). Similar to (fT8l ), terminations from group 3 
UEs to BS 2 cannot be considered to achieve the resultant connection state. Next, for at + b t = Mi, 

ct + dt < M 2 , b t < bt-i, c t < ct-i or a t + b t = Mi, c t + d t < M 2 , b t = b t -i, c t < c t -i or a t + b t = M u 



16 



a + dt < M 2 , b t < ht-i, c t = c t -i 

b t Ct-i 

(a t ,b t ,c t ,dt) _ (2) Q(b t ,ct,k,l)(2) Jb t ,ct,k,l) 



r 



k=0 1=0 



(a t ,k+b t+ i—b t ) (dt,l+c t+ i—ct) 
X (M.AO^LxA-x (M 2 ,A r 2 ) a d t _ 1 ,c t _i ' (21) 



The above case is an extension of the transition probability described in (I20I ). It considers the situation 
when full-cell occupancy occurs only at BS 1. Next, for at + b t = M\, ct + dt < M%, bt < bt-i, ct > Q-i 



r 



bt C t -! 

(a t ,bt,ct,d t ) _ (3)n(b t ,ct,k,l)(3){(b t ,c t ,k,l) ^ ^ n (a t ,k+b t+1 -b t ) ^ ^ n ( d *> 1 ) 



"if L t^t l f o tV"ti L tr)^ \u>t^~rvt+i—vtj yu't^; 

t-uct-i 56 t _ 1)Ct _ 1 (Mi,JV 1 )«o t _ 1 ,6 t _ 1 (Ma.JVa)^.,,^,- 0") 

fc=0 2=0 

The above case is an extension of the case described by (fT5l) . However, only full cell occupancy at BS 
1 is considered. Next, for at + h < M\, a + d t < Mi, h > bt-i, c t < c t -\ 

bt 

(a t ,bt,c t ,dt) _ (4)/ ) (6 t ,Ct,fc,Q)(4) t (b t ,ct,k,0) (a t ,k) (dt,c t -i—ct) rn-2\ 

'VxA-iA-iA-i ~~ 2^ V1A-1 ViA-i (^i,JV 1 ) Q! at _ 1 ,6 t _ 1 (M 2 ,iV 2 )arf t _ 1|Ct _ 1 • 
fc=0 



The above case is an extension of (fl6l) . In the case of (1231) full-cell occupancy does not occur in any of 
the cells, therefore the respective summation terms from (TT6l ) accounting for the edge case are removed. 
Also, there is only one summation because as in ( fT8l ) and d20b terminations from group 3 UEs to BS 2 
cannot be considered to achieve the desired end state. Lastly, for at + bt = M\, ct + d t < Mi, b t > bt-i, 
c t < C t -l 

bt-i c t 

(a t ,bt,c t ,dt) _ ST~^ ST~^ U) Q {bt,c t ,k,l)(A) Abt,c t ,k,l) (a t ,k) (dt,l+Ct-i—c t ) 

ViA-iA-iA-i-LL Via-. ?6 t _ 1 ,c t _ 1 (M 1 ,iV 1 )a 0t _ 1 ,6 t _ 1 (M 3 ,JV 2 )a (it _ ljCt _ 1 • (24) 
k=0 1=0 

The final case is an extension of the case described by (l23l . However, it considers the case when full 
cell occupancy occurs only at BS 1. 



E. Performance Metrics 

Given the complete description of the system, we are able to derive important metrics that would 
describe the efficiency of the load balancing process involving connection admission. While there are 
many performance metrics that can be extracted given the above framework, we focus on three primary 
metrics: (i) the blocking probability, which describes the probability that at least one UE which requests a 
connection from a particular group will be denied access to a channel, (ii) the channel utilization, which 
expresses the fraction of the available channels are being used, and (iii) the collision probability on a 
control channel, which provides the probability that at least one requesting connection will be lost due 
to a collision with another UE. 
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1 ) Blocking Probability: As used here, blocking occurs when at least one UE requests a new connec- 
tion, but cannot be admitted to any BS due to lack of available channels. Since each group has access to 
a different number of channels and can follow a different connection strategy, it is necessary to define 
separate blocking probability metrics for UEs in groups 1 and 2, as contrasted with UEs in group 3. For 
groups 1 and 2, the blocking probability is defined according to 

x N z y 

^=EE E E^,c, d T^S^ +k \ (25) 

a,b,c,d i=0 j=M v ~x-y+l k=0 

where for z = 1 x = a, y = b, v = 1 and for z = 3 x = d, y = c, v = 2. For group 3 UEs, the blocking 
probability is given as 

a Ni d N 3 b c N 2 
a,b,c,d ia.=0 j a =0 id=0 jd=0 i=0 j=0 k=cj> 

where <fi is defined separately for K < M2 and K = M2. For K = M2 <ft = M% + M2 — a — b — c— d — 
9a- 9d + i a + id + 1, where 



9a 



Mi - a + i a -b + i, j a > M\ — a + i a — b + i, 

(27) 

j a , otherwise, 



and gd is defined as (|27T ) replacing M\ with M2, a with c, i a with id, b with d and i with j. For if < M2 
<j) = M\ — g a — b — j a + i a + i + g c + 1, where g a is defined as (IITI) and 

max{0,sr e }, g e < K 

.'/<■ "U-c + j, c - j < if and ^ e > K, (28) 

0, otherwise, 
where g e = M 2 - c - d + i d + j - j d . 

We briefly explain the above equations. The derivation of the blocking probability for group 3 UEs is 
more complicated than for those UEs in groups 1 and 2 because this group can access channels from 
both cells. Therefore, the number of generations for group 3 UEs that leads to blocking has to account for 
the terminations within the same group, and also for the possible changes in the number of connections 
of UEs in groups 1 and 2. 

The blocking for group 3 UEs can be analyzed in two separate cases. The first case accounts for the 
number of generations needed to occupy all the channels in cell 1 (K = M2), while the second case 
accounts for the number of generations needed to occupy all available channels in cell 2 (K < M2). The 
first case is simpler to analyze because group 3 UEs need only generate as many connections as there 
are available resources on cell 1. 
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When K < M2, group 3 UEs can only access a maximum of K channels on cell 2. Therefore, extra 
conditions are added for the situation in which group 3 UEs are blocked when exceeding K connections 
in cell 2. If there are less than K available free channels after terminations of group 3 UEs connected to 
cell 2, the function max{0, g e } represents the number of connections to cause blocking by generating the 
exact number of connections to occupy all available channels. The max{0,<7 e } function is used to lower 
bound the necessary number of connections for blocking. This is because the number of connections on 
cell 2, in general, is not restricted to K and can thus have more than K current occupancies resulting 
in a possibly negative value for g e . On the other hand, if there are more than K available channels after 
terminations, exactly K channels are used by group 3 UEs to cause blocking. 

2) Channel Utilization and Total System Throughput: The overall channel utilization is 

U = E (« + 6 + c + ( 29 ) 

a,b,c,d 

which refers to the fraction of the collective capacity that has been used by all UEs in all groups. The 
average total system throughput is obtained by multiplying (l29l ) by H. 

3) Collision Probability: Because the model uses a random access channel for connection requests, 
it is necessary to compute the collision probability of the system. The collision probability for UEs in 
group x is: 

Nx-a K /AT \ 

d [x) = E E EtfMV^fV > fe d -P) N -"- k w%, oo) 

a,b,c,d k=0 j=0 ^ ' 

where r\ = {a, d} for group x = {1, 2}, respectively, and r\ = b + c for group x = 3. 

IV. Results 

Since our model incorporates a very large number of parameters, in the interest of clarity and brevity 
we focus our study on certain scenarios that are the most important in the context of our model. First, we 
present results that demonstrate the impact of a varying channel quality on the load balancing efficiency. 
Second, we examine the influence of the random access phase on load balancing efficiency. Lastly, we 
provide insight on the optimal channel sharing policy between BS 1 and BS 2. 

To confirm the correctness of the analytical model, we created a simulation environment for verifying 
the analytical results. The results in Sections IIV-AI and IIV-BI are obtained using both the analytical 
and simulation approaches to confirm correctness, while those in Section IIV-CI are obtained using only 
simulation. 
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Fig. 2. The channel utilization as a function of traffic intensity p. Two extremes of shared channels, i.e. K — (no load 
balancing) and K = M2 = 3 (all of cell 2's channels used in load balancing) are shown. The figure shows good agreement 
between the results from the analytical model and from simulation. 

A. Impact of Channel Quality on Load Balancing Process 

In this simulation, we model, among others, the call admission, termination, and load balancing 
processes exactly as described in our system model. As an example, we consider a scenario in which 
two identical cells are positioned such that they form a small TTR. For simplicity, we assume that 
N x = L x = 6, where x G {1,2,3} and Mi = M2 = K = 3. This particular analysis repre- 

(i) 

sents the effect of an increasing w% [ on the overall system-wide channel utilization, while setting 
Wi \ = = = 0.806 for all i G {u, d}, assuming reciprocal uplink and downlink conditions. This 
is equivalent to varying c^i from a location that is out of range to being right next to BS 1, and setting 

(i) 

dl,i = ds t 2 = ^2,2 = 30 m. We use the pathloss model with the following parameters: 7^ = — 85dBm, 
P t (i) = 30dBm, WW = -31.54 dB, a* = 3.65 dB, djj = lm and lastly 5 = 4.76. Furthermore, 
we assume an average channel capacity of H = 250 Kbps per channel with an average packet size 
r p = lkB and t s = 8 ms slot length. This yields an average packet length of l/q = 31.25 frames, or 
equivalently 250 ms. The channel throughput represents a typical value used in radio access network 
planning calculations [271 Table 8.17]. The packet size represents a realistic packet length sent over the 
Internet |[28l . where packets are distributed between a minimum value of 40 B (Transport Control Protocol 
acknowledgement packet) and a Maximum Transmission Unit, which for IPv6 equals 1.268 kB, for IEEE 
802.3 equals 1.492kB, for Ethernet II equals 1.5 kB, and for IEEE 802.11 equals 2.272kB. 

In order to determine the best performance of the load-balancing scheme, represented by the proposed 
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Fig. 3. The channel utilization, represented as a percentage on the vertical axis, as a function of the channel quality between 
group 3 UEs and BS 1, for two extreme values of K. As improves, more group 3 UEs generate successful connections 
to cell 1 resulting in more UEs that connect to BS 1 and consequently are offloaded onto cell 2, resulting in an overall increase 
in channel utilization. 

DR/< algorithm, for a varying channel quality we determine the level of traffic intensity which results 
in maximum channel utilization. Fig. |2] expresses the channel utilization as a function of increasing 
traffic intensity pi = P2 = P3 = P for two extreme K values, i.e. K = (when no load-balancing is 
used) and K = M2 = 3 (all of cell 2's channels may be used for load balancing). As expected, the 
channel utilization increases with more traffic intensity because an increase in p results in more frequent 
connection requests from UEs in all groups leading to a higher probability of successful connections. 
The increase in channel utilization tails off as the system reaches saturation, i.e. close to 100% channel 
utilization. Similarly, an increase in K results in a higher channel utilization as more UEs from group 3, 
that are blocked from cell 1, are offloaded onto cell 2 where they have access to an additional K channels. 
We observe that there is a decreasing rate of gain in channel utilization with an increase in the number 
of shared channels K. As Fig. [2] shows, a traffic intensity of p = 0.4 results in the most improvement 
in channel utilization between the two extremes of K = 0, and K = M2 = 3. With the knowledge of 
decreasing gains in channel utilization with increasing K, there may exist an intermediate value of K 
that not only leads to an improvement in total channel utilization, but also maximizes improvement with 
respect to the overall UE experience. This value of K is explored in Section ITV-CI In the current section 
we continue our investigation using p = 0.4 and explore the impacts of channel quality on performance. 



» K=0 (Sim.) 

. K=0(An.) 

o K=3 (Sim.) 
K=3 (An.) 
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Fig. [3] illustrates an increase in channel utilization with an increase in for two extreme values of 

(i) 

K. Increasing w\ \ results in group 3 UEs having more successful requests for receiving downlink 
transmissions because the average channel quality, in which requests are granted for group 3 UEs, 
improves. Therefore, ignoring the channel effects by assuming perfect channel conditions (also done 

(i) 

by setting [ = 1 in our model) in the analysis of load-balancing schemes, even for one particular 
group of UEs, produces a non-trivial difference in the channel utilization and leads to an exaggerated 
improvement in performance due to load balancing. By selecting a reasonable scheme to determine the 
channel quality, as presented in Section IIII-Al we are able to provide a more realistic evaluation of the 
improvements of load balancing. Note that the average channel utilization significantly increases as more 

(i) 

channels can be borrowed from BS 2. When \ increases, the difference in channel utilization between 
cases K = and K = 3 becomes more profound. This proves that with the low channel quality system- 
wide improvement from load balancing might not be as significant as in the case of perfect channel 
conditions. 

In Fig. g] we examine and B^ 2 \ as a function of for K = {0,1,2,3}. The 

scenario used in this result is identical to the one used previously. By comparing an increasing value 
of K in Fig. |4(a)j -Fig. |4(d)| we observe the impact of an increasing number of shared channels on the 
considered performance benchmarks. The first interesting observation is that irrespective of the value 
of K the blocking probability for group 1 UEs, is relatively stable. This means that the quality 

of service requirements for UEs not involved in load balancing will be met, even with load balancing 
enabled. Second, as K increases, blocking probability for group 3 UEs, B^\ significantly decreases, 
which proves the effectiveness of load balancing in this context. Furthermore, the collision probability for 
UEs in group 3, D^ 3 \ also reduces because with more shared channels there are fewer unconnected UEs 
to request connections. Note, however, that the difference in collision probability for an increasing K is 
not as significant as observed for the blocking probability. Finally, increasing K only slightly increases 
blocking probability B^> because these UEs have priority in connecting to any of cell 2's free channels. 

Focusing on Fig. |4(d)| only, where load balancing is enabled, we note that as increases, all curves 
experience an increase. This can be explained as follows: with an increase in w%\, more group 3 UEs are 
able to generate successful connections to BS 1 resulting in an increase in the contention for sub slots 
during admission control, and hence an increase in D^ 3 \ Also, there is an accompanied increase in 

^(3) 

because as more UEs generate successful connections, an increasing number of UEs contend for free 
channels on both cell 1 (where load balancing does not occur) and cell 2 (where load balancing occurs). 
Consequently, this also results in an increase in B^ and B^ 2 \ Although these trends are obvious, the 
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Fig. 4. Illustration of the relationship of the blocking probability for group 3 UEs, B^ 3 ', the collision probability for group 3 
UEs, D^ 3 \ the blocking probability for group 1 UEs, and the blocking probability for group 2 UEs, B^ 2 \ as a function 

of the channel quality between group 3 UEs and BS 1, Wg \ and (a) if = 0, (b) if = 1, (c) if = 2, (d) if = 3. We observe 
that the primary factor in the degradation of performance for group 3 UEs in this particular scenario is B^ 3 ' as compared to 
and this difference increases as improves. Once again, this figure shows good agreement between the results from 
the analytical model and from simulation. 



exact degradation in UE experience for each group is not. For example, in this specific scenario, Fig. [4] 
illustrates that is always the primary factor in the degradation of the group 3 UE experience as 
compared to D^ 3 \ This knowledge is significant as the network operator can determine whether an 
increase in K, or an increase L3 will be more beneficial to group 3 UEs. Observe that Fig. [3] and Fig. [4] 
show an extremely good match between the analytical result and simulation. 
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B. Impact of Random Access Phase on Load Balancing Process 

In this section we present results on the effect of random access phase on the performance of load 
balancing. The results are presented in Fig. [5] All network parameters are set identically to the network 
considered in Section HV-A1 except for L x = 3, where x G {1,2}. 

We begin by investigating the impact of different UE distributions on the performance of load balancing. 
We perform three experiments and denote each experiment as a specific case. In the first case we set the 
number of UEs, such that more UEs are distributed in groups 1 and 2, than in group 3, i.e. JVi = N2 = 6, 
N3 = 4. In the second case we set the number of UEs equal in each group, i.e. Ni = N% = N% = 6. 
And finally, in the third case we set the number of UEs in group 3 larger than in the other two groups, 
i.e. Ni = N2 = 6, A r 3 = 8. The metric that is studied in the three cases described above is the total 
network-wide blocking probability, calculated as | Yli=i as a function of channel access probability 
Pi = p for i G {1, 2, 3}. This metric is used in order to give a simple overall indication of the blocking 
suffered by UEs in all groups. Results are presented in Fig. |5(a)| 

The most interesting observation from Fig. |5(a)| is that with an increase in the number of UEs in the 
TTR, the total blocking probability becomes smaller for moderate values of p. Surprisingly, the blocking 
probability starts to drop sharply as the value of p continues to increase. This phenomenon occurs because 
as p increases, so do collisions on the random access channel, which in-turn limits the blocking probability 
because there are fewer UEs that successfully access available channels. It has to be kept in mind that 
for each case presented in Fig. |5(a)| the length of the random access phase remains the same. What is 
important to note is that for moderate values of p, the difference between blocking probabilities for each 
case is small, i.e. less than 5% (please compare values of blocking probability for each case in the range 
of p G (0, 0.6)). However, as p becomes very large, the curves with a higher number of group 3 UEs drop 
off faster because they experience a substantial increase in the number of collisions. Therefore, a certain 
value of p can maximize the channel utilization achieved through load balancing and also maintain the 
blocking probability at approximately the same level (given negligible changes in UE distribution). 

We now move our focus to the impact of random access phase length L x on the performance of load 



balancing. The results are presented in Fig. |5(b)| The set of parameters remain the same as in the earlier 
experiment in this section, however p x = p = 0.4. As an example, three network metrics are evaluated as 
a function of number of slots in the random access phase, L x : (i) total channel utilization in both cells, 
U, (ii) collision probability at group 3, D^ 3 \ and (iii) blocking probability at group 3, B^ 3 \ For clarity, 
the number of slots is set equal among each group of UEs. 
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Fig. 5. Impact of random access phase on load balancing process: (a) total network-wide blocking probability as a function of 
access probability p; (b) impact of random access phase length L x on the performance metrics of the considered system. We 
observe that the blocking probability is not a monotonous function of p and there is an extremum beyond which blocking starts 
to drop-off. As a result, the network metrics can be easily adapted by the network operator by dynamically selecting number of 
random access slots. 



Obviously, as the number of random access slots increase the collision probability decreases for group 
3 UEs, and the overall channel utilization increases. However, as the collision probability decreases the 
blocking probability, within the same group of UEs, becomes larger. This is in line with the results 
presented in Fig. |5(a)| Recall, that as more UEs gain access to the BS, the probability that channels 
become unavailable increases. The results shown in Fig. |5(b)| further demonstrate the fundamental tradeoff 
between the delay caused by random access and overall network utilization. Apart from this tradeoff, 
we demonstrate that the network operator has a powerful tool, i.e. random access phase length, through 
which network metrics can be easily regulated. It is obvious that the operator has no control over the 
channel access probability, p, of individual UEs. However, the operator is able to set a higher value of 
L x to the ASC of interest in order to maintain an expected access delay for each UE against a required 
blocking probability. 

C. Impact of Varying Shared Channel Pool K on Load balancing Efficiency of DRk 

We consider a macrocell scenario in which the distribution of UEs in groups 1, 2 and 3 follow the 
relationship given by N\ = N2 < N3 and L x = N x , where x E {1,2,3}. Let N\ = N2 = L\ = 
L2 = 25 and N3 = L3 = 40. We consider a symmetric system where each cell has 10 channels, i.e. 
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Fig. 6. The percentage improvement of channel utilization (left vertical axis) and the percentage improvement of channel 
utilization per shared channel (right vertical axis) as a function of shared channels, K. As K increases, there is an obvious 
improvement in channel utilization, however, there are decreasing gains experienced per additional K. 

M\ = M2 = 10, and the distances between each group of UEs and their respective serving BSs are 
identical, i.e. d% t \ = ^31 = cfe^ = c?2,2 = 220 m. Once again, we assume an average channel capacity 
of H = 250 Kbps per channel with an average packet size r p = lkB. A frame duration duration of 
t s = 1ms is used. We assume a simplified path loss model with identical parameters as in Section IIV-AI 
except with 5 = 3, which is more appropriate for outdoor channel conditions. 

In Fig. [6] we explore the system- wide improvement in channel utilization (represented as a percentage 
on the vertical axis) as a function of the number of shared channels K (represented on the horizontal 
axis). Note that for the remainder of our study we fix traffic intensity for all groups to p = 0.2 in 
order to determine the maximum gain in channel utilization for every value of K. The line with circle 
markers represents the total percentage improvement experienced as a function of K, while the line 
with triangle markers represents the improvement experienced in channel utilization per shared channel. 
As K increases, there is a decreasing improvement in channel utilization per additional channel, which 
indicates that the extra cost of sharing more channels for load balancing may outweigh the added benefit 
of serving a greater number of UEs. 

Although there is an overall improvement in the system-wide channel utilization, the exact effect of the 
load-balancing scheme on the UE experience is unknown. Obviously, decreases with an increase in 
K because group 3 has access to more channels. In contrast, B^> increases because more UEs in group 
3 access channels belonging to cell 2, which are of course also accessible to UEs in group 2. Although 
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Fraction of channels in BS 2 (K) 

Fig. 7. The total change in blocking probability, B^ x \ and blocking probability per shared channel, B^ x ' /K, where x £ {2, 3}. 
In this scenario, we observe that the decrease in is always greater than the corresponding increase in B (2) for all K, 
suggesting that there is an overall improvement in the UE experience. Furthermore, it is seen that at K = 6 we have the greatest 
difference between the decrease in B^ /K and the corresponding increase in B^ /K suggesting that this is the optimal number 
of shared channels to use in order to gain the best UE experience per K. 

these general trends are obvious, the exact relationship between the amount of decrease in versus 
the amount of increase in B^ is unknown. In Fig. [7] we examine this relationship in more detail, where 
the decrease in B^ (solid line with triangle markers) is plotted with the consequent increase in 
(solid line with circles) as a percentage on the vertical axis with increasing K on the horizontal axis. We 
observe that in this particular scenario the total decrease in B^ is always more than the total increase 
in B( 2 \ suggesting that the overall UE experience improves with the proposed load-balancing scheme. 
This reaffirms the increase in channel utilization with an increase in K for DR^, which is previously 
observed in Fig. [6] 

The total improvement in overall UE blocking probability demonstrates the effectiveness of the load 
balancing scheme. However, from a network operator viewpoint, knowledge of the changes in UE 
experience per additional shared channel is also very important. Fig. [7J examines the effect of an increase 
in K on both the decrease in B^/K (dashed line with triangles), and the consequent increase in 
B^/K (dashed line with circles). We observe, that for this particular scenario, a decrease in B^/K is 
always more than the increase in B^/K, which suggests that the UEs in group 3 experience more of 
an improvement in performance than the performance degradation experienced by UEs in group 2 per 
additional shared channel. This allows direct evaluation of the effectiveness of the load-balancing scheme 
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on the overall UE experience per additional shared channel. In Fig. |7J we note that B^/K reaches a 
maximum at an intermediate value of K, i.e. K = 3, because the system reaches a balance between 
the number of UEs requesting connections and those that are already connected through load balancing. 
Our model allows for the direct observation of this system state because of the combined modeling of a 
finite number of UEs together with a detailed call admission process. With the use of Fig. |7J we are able 
to determine the best K to improve the overall UE experience on a per shared channel basis, and then 
find the corresponding improvement in the overall channel utilization using Fig. [6l For this particular 
scenario, the maximum difference between the increase in B^/K and decrease in B^/K occurs at 
K = 6 and corresponds to an overall improvement in channel utilization of 12.6%. 

In summary, we can construct the following optimization function. Given fundamental descriptors of 
the network considered in Section JIIJ i.e. 2Vj, Mi, L x , pi, Wx} y , r p , and t s , the network operator should 
find 

argmaxU/K subject to ViB^ < 6 (i) ,Z? W < d (i \ (31) 

K 

where bi and di are the required maximum blocking and collision probabilities, respectively, for group 
i. The developed analytical model provided in Section HIT] allows solving the optimization function < f3TT >, 
since each metric, U, B^ 1 ' and is given in closed form. The optimization formula allows obtaining 
the value of K in order to obtain the maximum utilization per shared channel, such that all considered 
quality of service metrics required by the operator, bi and di, are met. Note that finding the optimal 
solution to (f3TT > is beyond the scope of this paper. 

V. Conclusions 

We have presented a new analytical model to assess the load balancing process in cellular networks. 
The model differs in many respects from previous work on load balancing analysis in several important 
ways. In particular, it incorporates a detailed call admission procedure, addresses traffic prioritization in 
the connection admission process, and also allows exploration of the impact of channel quality on load 
balancing efficiency. We have presented a variety of results derived from this framework, and in particular 
explored the tradeoffs in terms of channel utilization, blocking probability, and collision probability when 
traffic is transferred from a highly congested cell to a less-loaded neighboring cell. 
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